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Abstract. Analyticity and unitarity techniques are employed to obtain bounds on the shape 
parameters of the scalar and vector form factors of semileptonic K13 decays. For this purpose 
f> I we use vector and scalar correlators evaluated in pQCD, a low energy theorem for scalar form 

D . factor, lattice results for the ratio of kaon and pion decay constants, chiral perturbation theory 

' calculations for the scalar form factor at the Callan-Treiman point and experimental information 

on the phase and modulus of Ktt form factors up to an energy ti n = lGeV 2 . We further derive 
regions on the real axis and in the complex-energy plane where the form factors cannot have 
zeros. 



1. Introduction 

K13 decays are important for determining the matrix element V us , for recent reviews, see [HE] 
of the Cabibbo-Kobayashi-Maskawa (CKM) which in turn is crucial for testing the unitarity of 
the CKM matrix. Information on the experimental sector is rich, in particular for the Kg^ decay 
rates, which were measured by a number of experiments BNL-E865 [3], KLOE [3], ISTRA+ [5], 
KTEV [B] and more recently a new analysis by NA48 [7]. Recent lattice studies have also been 
carried out, see refe.® El EDI EE]- 

The decay of a kaon to a pion, a charged lepton and a neutrino is described by the matrix 
element 

(n°(p')\s^u\K + (p)) = + (p-pV_(t)], (1) 

where /+(£) is the vector form factor and the combination 

M) = U® + Mf^f-W (2) 

is known as the scalar form factor. The expansion at t = 

/*(*) = /+(°) ( l + ^if2 + + • • • ) , (3) 



defines the slope A^ and the curvature A^ parameters where k = denotes the scalar and k = + 
denotes the vector form factor. The precise determination of the element \V US \ depends on how 
accurate the parametrization of the form factor is. To improve the precision and to provide 
bounds on the shape parameters of the form factors, we use inputs coming from certain low 
energy theorems, perturbative QCD, lattice computations and chiral perturbation theory. Our 
techniques allow us to incorporate the phase and modulus information of the form factors. We 
also apply the technique to find regions on the real axis and in the complex t-plane where zeros 
are excluded. The knowledge of zeros is of interest, for instance, for the dispersive methods 
(Omnes-type representations) and for testing specific models of the form factors. For more 
details, see refs. [TJ CE3J CES CE5] . 



2. Formalism 



The formalism described in [121 [TBI [14"] and in the contribution to these Proceedings [16] , 
exploits an integral of the type 

poo 

/ dt P+fi (t)\f +fi (t)\ 2 <I +fi , (4) 
Jt+ 

along the unitarity cut, whose upper bound is known from a dispersion relation, satisfied by a 
certain QCD correlator. For the scalar form factor this reads 

d r 2tt , I f°° 7 tlmlWt) 

Imllo(t) > p l/owl > I 6 ) 

withi± = (M K ±M 7T ) 2 . Analo gues expression, involving a suitable correlator denoted by XiiQ 2 )-, 
can be written down for the vector form factor. We can now use the conformal map t — > z(t) 

that maps the cut t-plane onto the unit disc \z\ < 1 in the z = zit) plane, with t+ mapped onto 
z = 1, the point at infinity to z = — 1 and the origin to z = 0. Using this map, we cast eqn.Q 
into a canonical form, incorporate phase and modulus information as well as the Callan-Treiman 
relations and finally employ a determinant for obtaining bounds on the shape parameters and 
for finding regions of excluded zeros in the complex t-plane. 



3. Inputs 

The essential inputs of our formalism, the vector Xi{Q 2 ) an d the scalar Xo(Q 2 ) correlators, 
can be calculated in perturbative QCD up to order for Q 2 » A.q CD [T71 [JS]. We get 
Xi(2GeV) = (343.8 ± 51.6) x 10~ 5 GeV~ 2 and xo(2GeV) = (253 ± 68) x 10" 6 , see also [T2]. 
An improvement can be achieved when we implement theoretical and experimental information 
into the formalism of unitarity bounds. The first improvement comes when we use the the value 
of vector form factor at zero momentum transfer. Recent determinations from the lattice give 
/+(0) = 0.964(5) |19| . We can also use two low energy theorems, namely soft pion and soft 
kaon theorems, for the improvement of the bounds on the slope and curvature parameters in the 



scalar case. The soft pion theorem relates the value of scalar form factor at first Callan-Treiman 
piont A Kir = - M 2 to the ratio Fk/F^ of the decay constants [201 12T]: 

f Q (A Kn ) = F K /F n + Act- (8) 

Recent lattice evaluations with Nf = 2 + 1 flavors of sea quarks give Fk/F^ = 1.193 ± 0.006 
[H [9]. In the isospin limit, Act = —3.1 x 10~ 3 to one loop [22] and Act — to two-loops in 
chiral perturbation theory [23 1 124 1 f25] . 

At Ak%{= —Aktt), a soft-kaon result [26] relates the value of the scalar form factor to F^/Fk 

fo(-A Kn ) = FJF K + Act- (9) 

A calculation in ChPT to one-loop in the isospin limit [22J gives Act = 0.03, but the higher 
order ChPT corrections are expected to be larger in this case. As discussed in [J2], due to the 
poor knowledge of Act 5 the low-energy theorem eqn.Q is not useful for further constraining 
the shape of the form factors at low energies. On the other hand, we obtain from the same 
machinery, the stringent bound on the quantity Act which is —0.046 < Act < 0.014. 

Further improvement of the bounds can be achieved if the phase of the form factor along 
the elastic part of the unitarity cut is known from an independent source. In our calculations 
we use below t{ n the phases from [271 [28] for the scalar form factor, and from [291 [30] for 
the vector form factor. Above t; n we take 6(t) Lipschitz continuous, i.e., a smooth function 
approaching tt at high energies. The results are independent of the choice of the phase for 
t > t- m . We can further improve the bounds if the modulus of the form factor is known along 
the unitarity cut, t < t; n : we can shift the branch point from t± to ti n by subtracting the low 
energy integral from the integral Eq. (jl]). In order to estimate the low-energy integral, which 
is the value of the integral contribution from t + to t{ n , see for expressions ref.[13], we use the 
Breit-Wigner parameterizations of and |/o(i)| in terms of the resonances given by the 

Belle Collaboration [31] for fitting the rate of r — > Knu decay. The above leads to the value 
31.4 x 10~ 5 GeV~ 2 for the vector form factor and 60.9 x 10~ 6 for the scalar form factor. By 
combining with the values i+,0, we obtain the new upper bound on the integral Eq. (j4|) from 
t in to oo, I' + = (312 ± 69) x 10" 5 GeV- 2 and I' Q = (192 ± 90) x 10" 6 . 

4. Results 

In Figs. [1] and El the constraints for the scalar form factor are represented together with 
experimental information from various experiments. As shown in Fig. [H the slope Aq of the 
scalar form factor, predicted by NA48 (2007) is not consistent with our predictions (yellow band) 
which are obtained by taking into account the phase, modulus as well as the CT constraint. 
Nevertheless our predicted range for the slope is well-respected by the recent 2011 analysis by 
NA48 [7J . 

The value of Aq for this new determination by the NA48 reads A = (15.6 ± 1.2 ± 0.9) x 10~ 3 . 
On the theoretical side, the prediction of ChPT to two loops gives Aq = (13.9~?;| ± 0.4) x 10" 3 
and Ag = (8.0+J;^) x 10" 4 which are consistent with our results within errors as shown in Fig. [2J 
For the central value of the slope A given above, the range of A ' is (8.24 x 10 -4 , 8.42 x 10 -4 ). The 
corresponding theoretical predictions are A = (16.00 ± 1.00) x 10~ 3 , Aq = (6.34 ± 0.38) x 10~ 4 
obtained from dispersion relations. 

Comparison of the experimental results with our constraints for the vector form factor is 
shown in Fig. [3l We find that except for the results from NA48 (2007) and KLOE, which have 
curvatures slightly larger than the allowed values, the experimental data satisfy our constraints. 
The new results from NA48 [7] provide a curvature which overlaps with our constraints while the 
slope lies completely within our domain as seen in the figure. We also note that the theoretical 
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Figure 1. The allowed range for the slope of the scalar form factor, when we include phase, 
modulus and the CT constraint (yellow band). The grey band shows the range without the CT 
constraint. 
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Figure 2. Allowed domain for the slope and curvature of the scalar form factor, using the 
normalization /+(0) = 0.962, the value fo(AKn) = 1-193, and phase and modulus information 
up to t in = (lGeV) 2 . 



predictions A' + = (24.9 ± 1.3) x 10" 3 , A" = (1.6 ± 0.5) x 10~ 3 obtained from ChPT to two 
loops, and A' + = (26.05^;^) x 10" 3 , A'j = (1.29_tg;{£) x 10" 3 , and A' + = (25.49 ± 0.31) x 10~ 3 , 
A" = (1.22 ±0.14) x 10~ 3 obtained from dispersion relations are consistent with the constraint. 
For more results, see |12| . 

We can also extend our technique to derive regions in the complex plane where the form 
factors can not have zeros. For the Ktt form factors, the influence of possible zeros in the 
context of Omnes dispersive representations has been analyzed in |30j. The absence of zeros is 
assumed in the recent analysis of KTeV data reported in [BJ. In Fig. H]we show the region where 
zeros of the scalar form factors are excluded. The formalism rules out zeros in the physical 
region of the kaon semileptonic decay. In the case of complex zeros, we have obtained a rather 
large region where they cannot be present. In the case of the vector form factor, the analysis of 
|30j using data from r decay concludes that complex zeros cannot be excluded, due to the lack 





Figure 4. Domain without zeros for the scalar 
form factor: the small domain is obtained 
without including phase and modulus in the 
elastic region, the bigger one using phase, 
modulus and CT constraint. 



Figure 5. Domain without zeros for the vector 
form factor: the small domain is obtained 
without including phase and modulus in the 
elastic region, the bigger one using phase and 
modulus. 



of information on the phase of the form factor in the inelastic region. However our formalism 
is independent of phase information in the inelastic region and leads without any assumptions 
to a rather large domain where complex zeros are excluded. For the vector form factor, fig. [5] 
shows the region where complex zeros are excluded. For more results, see |12j . 



5. Conclusion 



We have derived stringent constraints on the shape parameters of the form factors of 
decay which is the best source for the extraction of CKM matrix element V us . The results 
are promising and stringent especially in the case of the scalar form factor. The most recent 
results from NA48 [7j is consistent with our prediction for the slope of scalar form factor and 
restricts the range of the slope to ~ 0.01 — 0.02. We have also excluded zeros in a rather 
large domain at low energies both for the scalar and vector form factor. The Callan-Treiman 
input provides an additional constraint in the case of the scalar form factor and as a result 
excludes a larger domain of the energy plane where zeros can exist. Thus, this work represents a 
powerful application of the theory of unitarity bounds, which relies not so much on experimental 



information, but on theoretical inputs from perturbative QCD, low energy theorems and lattice 
calculations. It provides a powerful consistency check on determinations of shape parameters 
from phenomenology and experimental analyses. 
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